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Abstract 



We calculate the as corrections to the form factors which parameterize the hadronic 
tensor relevant for inclusive semileptonic B XrVj- and A;, —* XtVj- decays. We apply 
our results to the double differential decay rates for the decays of B mesons and polarized 
Af, baryons to polarized r leptons, presenting them in terms of one-dimensional integrals. 
Formulas appropriate for insertion into a Monte-Carlo simulation are presented. 
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1. Introduction 



There has been much recent interest in inclusive semileptonic B decays as a means of 
measuring the CKM parameters V^b and Vcb- These decays are very difficult to calculate 
due to quark confinement and the breakdown of perturbation theory. However it is possible 
to improve the behavior of the QCD perturbation series by considering a "smeared" decay 
rate, in which one integrates in a certain way over the lepton momenta. The infrared 
singularities responsible for the breakdown of pertubation theory are tamed by integrating 
appropriately across perturbative thresholds. The idea of smearing to improve the behavior 
of the QCD pertubation series in the physical region has been used for example to examine 
the inclusive process e+e" hadrons ^ For this process one computes averages of the 
QCD contribution to the vacuum polarization. This smearing method can also be used 
to describe inclusive semileptonic heavy meson decays B XlV [^. In this case there 
are additional non-perturbative complications, since one must compute the expectation 
value of quark currents in a S meson rather than the vacuum. This difficulty is controlled 
by means of a heavy quark expansion. In this paper, we compute perturbative QCD 
corrections to semi-inclusive semileptonic heavy meson decays involving tau leptons, in 
the framework of the above smearing method. We will consider the leading term in the 
heavy quark expansion, so that in practice no non-perturbative matrix elements will be 
needed. 

The semileptonic B meson decays are caused by the Weak Hamiltonian density 

AGf. 
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Hw = -V,b^q,l^PLhht.PLi^. (1-1) 



where Pl is the chiral projector ^(1 — 75). The index j denotes either a charm or an up 
type quark. The inclusive decay rate is given by 



dq^dE^dEr ~ 27r3 



2 n2 

L^^W^r (1.2) 



L^i, is an easily computed matrix element of leptonic weak currents 

L^, = Tr[{i)i + 77107^(^)7.^^1.], (1-3) 
and Wfj_jj is the hadronic matrix element 

W^, = (27r)3 Y^5\p^-q-p^) {B{v, s)\ \X) {X\ J/ \B{v, s)), (1.4) 



X 



with the weak quark current Uj'f'^PLb and q = pi + Pu the total lepton momentum. 
Using the optical theorem, W^'^ may be expressed as the discontinuity across the cut in 
the time ordered product of weak currents in a B meson: 

H^M^ = -i/^T'^^^ (1.5) 

TT 

where 

T^^ = j d'^xexpiiq-x) (S| T[J^\x) J^(0)] \B) . (1.6) 

The B quark matrix element may be computed using a heavy quark expansion in — 

In the approach of ref [|| , there are two steps to this expansion. The first consists 
of an operator product expansion in which one expands about an "on-shell" value for the 
heavy quark momentum: p = mtv + k where mi, is the mass parameter of the heavy quark 
effective lagrangian, v is the four velocity of the B meson, and /c is a small residual off-shell 
momentum. The Wilson coefficients are written as a pertubation series in ctg. The second 
step is a heavy quark expansion of the expectation values of the local operators in a B 
meson. The Wilson coefficients are singular at the perturbative threshold, at which the 
invariant mass of the hadronic decay products is rrij . However, by appropriately smearing 
the decay rate over lepton energies the operator product expansion can be made well 
behaved simultaneously with the pertubative QCD expansion. This will be discussed in 
more detail in section 2. The lowest order term in the expansion is simply the parton 

model result for the appropriately spin averaged decay of a 6 quark. 

It has been shown that the corrections vanish [||, which lends considerable 

credence to smeared parton model calculations. The order (— f^)^ corrections have been 
computed at tree level 0, and the resulting parametrization has been applied to the 
differential decay rate for B —>■ XtV with polarized and unpolarized tau leptons. To 
order as the parton model result for semileptonic decays involving an electron have been 
calculated For decays involving a massive lepton in the final state, parton model as 

corrections to the total decay rate of S ^ X^tV have been computed for vanishing 
final quark mass and unpolarized tau lepton. Corrections to the differential decay rate ^ 
for B mesons going into polarized or unpolarized tau leptons have also been computed, for 
massless or massive final quarks, in the rV rest frame [|] . 

In this paper we compute the matrix element of time ordered weak currents in "off- 
shell" b quark states to order as- We apply this to double differential rates for B meson 
and baryon semi-inclusive semileptonic decays involving either massive or massless final 
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quarks, and either polarized or unpolarized tau leptons. The off-sheU result for the time 
ordered product is useful should one ever wish to calculate the as corrections to the next 
to leading terms in the expansion. We leave this for future work. By computing 

the absorbtive part of the on shell value, we arrive at the parton model results for decays 
involving a tau, in terms of a one dimensional integral over the neutrino energy. This decay 
rate differs from that involving an electron due to the large r mass. After a short review 
of the operator product expansion and the smearing method in section 2, section 3 studies 
the kinematics of B decays, section 4 gives our result for the B meson triple diferential 
inclusive decay width, section 5 studies the analytic structure of the form factors, section 
6 deals with divergence cancellations, section 7 gives the formulas for polarized decays, 
and section 8 states our conclusions. Our results for B decays are summarized, in a form 
appropriate for inclusion in a Monte Carlo simulation, in Appendix E. 



2. The Smearing Method 

In this section we briefly review the procedure of Chay, Georgi, and Grinstein by 
which one arives at the "smeared" parton model as the leading term in a expansion. 
The initial step in this procedure is an operator product expansion of the time ordered 
product of quark weak currents in terms of local operators in the heavy quark effective 
theory. 

f d^xexp{iq-x}T (^J^{x)r\o)') = ^(q, t;)C>,,„(0). (2.1) 

The local operators Oi^y{x) are made up of light degrees of freedom (gluons and light 
quarks) and the b heavy quark operator h'^{x), related to the full QCD field by 

hlix) = ^^exp{imi,yx)b{x), (2.2) 

with V the four velocity of the B meson. The coefficient functions Ci{q,v) can then be 
obtained to finite order in ctg by considering the perturbative matrix elements of both 
sides of equation (|2.1|) between (unphysical) free quark and gluon states, where the b 
quark states carry momentum rribV + k. The quark mass used here is the mass parameter 
of the heavy quark effective lagrangian. All the dependence on the residual momentum k 
is contained in the matrix elements of local operators containing covariant derivatives of 
the heavy quark field. The resulting operator product expansion is thus an expansion in 
off- shell momentum k. 
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After the expansion has been computed, it is reinserted between B meson states. 
One might hope to calculate or at least paramaterize the hadronic tensor in this way. 
However the operator product expansion breaks down near the perturbative threshold. 
Since the perturbative analytic structure of T^*^ is contained in the coefficient functions, 
the coefficient functions must be singular near thresholds. To see this explicitly, consider 
the tree level expectation value of the current-current correlator in a 6 quark: 



{rribV — q + kY — + 



u — 7 -TT^ — : u. (2.3) 



The tree level coefficient functions for two quark operators are then obtained by ex- 
panding the above expression as a Taylor series in k. The tree level operator product 
expansion has been considered previously by several authors [|] and the reader is referred 
to them for details. The important point is that the expansion parameter is 

k 



(2.4) 



The denominator vanishes at threshold, so the operator product expansion breaks down in 
its vicinity. The discontinuity across the real axis in the complex v-q plane consists of an 
infinite series of derivatives of delta functions which appears to be useless at finite order: 

= 5^ c„5(-) {{m.v - qf - m]) , (2.5) 

n 

where the are matrix elements of the general form 

Cn = {B\'^D^hl\B). (2.6) 

The perturbative coefficient functions possess singularities in the wrong places: to arbitrary 
perturbative order in ctg the threshold singularity is that of the parton model of free quarks 
and gluons, (m^f — q)"^ = m?j, rather than a real hadronic threshold. 

These troubles may be avoided by computing Ti(v-q, q^) well away from the physical 
region. A large imaginary piece to v-q acts as an infrared cutoff, allowing one to compute 
the Wilson coefficients perturbatively 0]. Furthermore, the coefficients get progressively 
smaller for higher orders in the residual momentum k. Far from the physical region, one 
expects both the operator product expansion and QCD pertubation theory to be well 
behaved. Thus one can use them to compute integrals of VFi(v-q, g^) over real v-q. Due 
to equation ( p..5| ), such an integral is equivalent to an integral of Ti(v-q, q^) just above the 




wwwv 



Fig. 1 Analytic structure of T^*^ and integration contour in 
the complex v ■ q plane 



cut minus an integral just below. The integration contour can be deformed away from the 
real axis between the endpoints of the integration. At a distance mf, from threshold in 
the complex v-q plane the effective expansion parameter becomes As long as these 
endpoints are far away from both perturbative and real thresholds, the error due to the 
integration near the endpoints should be small. 

In this paper we compute the differential decay width 'd<fdE- - These three kine- 
matic variables may be varied as long as we treat the final state quark as off-shell. Pertuba- 
tion theory does not give a sensible result for this rate, but one can obtain sensible results 
for -j^-^, provided that and are such that the boundaries of the v-q integration 
are not near threshold. Integration over v-q = + -E'77 is performed by integrating the 
triple differential decay width over neutrino energies at fixed tau energy. 

In the parton model, the cut in the complex v-q plane is given by the B decay threshold 

m? + — m'j 

There exists another cut in v-q, coming from the threshold for weak processes involving 
two b quarks in the final state. This threshold is given in the parton model by 

,.^^ i2m,^^.f-,nl-,^ ^ (2.8) 

2mi, 

For the particlar case rrij = and = m^, this threshold coincides with the B decay 



threshold (|2.7|). In this case, near that point one can not get the operator product expansion 
to work by integrating past the B decay threshold, because there is no way to close 
the integration contour in the complex plane far from the region where the perturbative 
expansion fails. The cuts together with the integration contour are drawn in fig. 1. 
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We shall compute the B decay rate using only the lowest order term in the operator 
product expansion. This calculation is equivalent to the parton model calculation of the 
spin averaged decay rate of a free b quark, where the quark mass is equal to the mass 
parameter of the heavy quark effective theory. The reason for the equivalence is as follows. 
The leading term of the operator product expansion is determined from the expectation 
value of the current correlator in an on shell b quark. This expectation value is of the 
form Cj{y-q, q^, as)u{mbv)Tju{mbv), where Tj is some combination of Dirac matrices and 
u is an on shell spinor. The corresponding term in the operator product expansion is 
Cj{v-q,q'^,as)hlTjhl ^. This last statement is true to all orders in as due to the heavy 
quark spin symmetry: at leading order in ^^^^ any operator of the form hlTjhl is a 
component of a conserved current associated with heavy quark number or heavy quark 
spin. Furthermore the expectation value of such an operator in a 5 meson is exactly 
known by the very same symmetry arguments; the Isgur-Wise function at zero recoil is 
one. Thus the B meson expectation value averaged over polarizations is equal to the b 
quark expectation value averaged over spins: 

T"" = ^a(v-q, q^)Tr (^F,^ , (2.9) 

where we have normalized our heavy quark states to v^. 



3. Kinematics 

It is convenient to decompose the hadronic tensor Tf^^y into form factors, 

= -gi^-'Ti + v''v''T2 - ie'^'^'f^VaqisTs 

^q^q'^T^ + {q^v" + q''v^')T^ (3.1) 
+ {q''v'' -q''v'')TQ. 

Tq vanishes due to the time reversal invariance of the strong interactions: T^'^{v,q) = 
T^^iy'q') where v'^ — and q'^ = q^. W^y may be similarly decomposed. The decay 
rate for spin averaged S's is given by 



|2 ^2 



dq^dE^dEy 



27r3 



(g2 - ml) + W2 ( 2E,Er - 



mz 



+ {{Er - £;.)(g' - ml) - 2E,ml) 



(3.2) 



+ ml [W^i 



2 1 
.q — 



)+Wr,{2E,] 



Operators which contain gluons or light quark pairs are sub-leading. Thus, weak annihilation 
of the b quark and light antiquark is suppressed in the mi, 00 limit. 
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For ± helicity tau's we have: 



d^T± 



1 



dq^dErdE^ 2 dq^dErdE^ 

|2 

I l/„-f. 

± 



1 ± 



E^ 



E. 



El 



47^3 



-2M^i -W2- IWz (^^E^ - ml cos a^E, 

+ VF"4 (m^ + 2Er{\/El - ml co-^oi- Er 
+ {2{^El-ml cos a- 



where the angle between the tau and neutrino trajectories is given by 

ml-q^-V 2E^Er 



cos a 



2E,^El-ml 



(3.3) 



(3.4) 



Note that this expression reduces to eq. ( |3.2| ) when summed over tau hehcities. The 
leptonic variables and Ejj = v-q — E^- are constrained to the region given by 



m^< E^< ^ 



"(1 



m 



B 



^ ^ (Er + J El - ml){m\ - m\ - ^E^ms) + m^ms 



ms — Et — y^El — ml 



Er + - 



- ml 



< v-q < Min[E^ + - 



m| + q'^ 



m 



X 



(3.5) 

(3.6) 
(3.7) 



2{E^ + y/El - ml) - ■ ■ 2{Er - ^E^ - ml) ' 2mB 

ms is the B meson mass, and mx the mass of the lightest hadron containing the 
final state m or c quark . With the usual assumption ms > ttii, , the allowed region of 
integration for v-q includes the parton model theshold ( p.7|) . The upper limit for v-q shows 
a new feature of tau kinematics; for a massless final lepton the upper limit is simply the 
threshold value (|2.7|) . The lower and upper bounds in (|3.7|) correspond to parallel and 
anti-parallel r and Uj- trajectories. 

In the case of a decaying polarized A;, we can parameterize the spin-dependent part 
of T^'^ in terms of nine additional form factors, 



-s-q 



-g^^'Gi + v^v''G2 - le^^'^^f^v^qpCs + q'^q'G^ 



+ {q^v'' + q'-vn + {s^v" + s^v^) Ge + {s^q" + s^q^) Gj 



(3.8) 
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where s is the spin vector of the A;, baryon in the rest frame. Note that there are no 
terms of the form i{q^e'^°'^'^Vaqpsx — (// ly)) or i{v^e"°'^^Vaqf3S\ — (// ^ u)). They can 
be shown to vanish using the identity 

l/ l/ 

The differential decay rate in this case has additional terms that depend on the spin of the 
A,: 



1 



dq^dErdE^dcosO 2 dq^dErdE^ 



+ cos 6 



1 rf^rf'-^-) 

2 dq^dErdE^ 



{Ej, cos a + ^JE"^ — ml) 



+ 



87r3 



1 ± 



- G7 {2mlE^ cos a) - Gg (2E^{Er cos a - ^/E^ - ml 

- Gq (^{q^ - ml){E„ cos a - A/Ef^^ml" ) + 2m^E^ cos a 

± — ^ = ( G% (cos a) + G7 ( E'r cos a — ^/ E^ — m 
^JEl - ml^ ^ 

+ Gs (cos a) + Gg (cosct(£'iy + -E't)) 



(3.9) 

where 6 is the angle between the spin vector of the Aj, and the r three-momentum, and 
where in the second term d^V^^''^^''^ is identical to d^V±, except that all the Wis are 
replaced by the spin-dependent G^'s. 

Note that when we integrate the polarized cross section ( p.9|) over cos^ the term 
proportional to cos 6 vanishes, while the constant term gets multiplied by a factor of 2 and 
so reproduces eq. ( p.3|) . 



4. B meson decay 

To compute the time ordered product of eq. ( |1.5| ), we must evaluate the graphs in 
figs. 2,3,4 and 5. We do this by first separating the Dirac structure from the loop integrals, 
evaluating the loop integrals in terms of a function we call S, and reducing the Dirac 
algebra with the help of the heavy quark effective theory. We present analytic formulas for 
S and express the form factors Ti in terms of S and its derivatives. After that, it is a matter 
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of numerical integration to generate plots at arbitrary values of and E^, with kinematic 
cuts, or with additional refinements. Such refinements might include corrections due to 
the W propagator or rrij/mb corrections arising from expressing the width in terms 
of quark masses evaluated at one scale instead of on-shell 0. In addition, higher order 
corrections due to these graphs (such as the correction to the operator {B\ (i-D)^ 

h 

may be computed by substituting p rritv + k and expanding in Strong corrections 
to decays of polarized A[,s may also be found straightforwardly by evaluating the Dirac 
traces between spinors instead of averaging over spin. 

We do our calculation in minimal subtraction and Feynman gauge, with a gluon 
mass A to regulate infrared divergences. For a massless final state, there is considerable 
simplification if one uses dimensional regularization instead of a gluon mass fl^, but for 
B clV, the simplification is less obvious. For the charm case, the amplitudes may 
be expressed in terms of generalized hypergeometic functions, which can be numerically 
integrated over the charm momentum and analytically integrated over the lepton momenta 



11}] . However, if lepton spectra are desired, this trick for analytic integration over lepton 
momenta is not useful. 

We devide our calculation into several parts. The total form factor is expressed as the 
sum of the contributions of various diagrams 

rp rpO _|_ rpSelf rpyert _j_ rphox 

4.1. Heavy Quark Self Energy corrections 

The quark matrix element of the tree level graph of fig. 2 is 

Mp) ., \ 2^ M p)^ (4-1) 

where p is the momentum of the external b quark, q is the W momentum, and Ub{p) is 
the b quark spinor. For the parton model result, we will eventually take p = mi,v, with 
v^ = l. 

Wave-function renormalization of the external heavy quark lines modify the tree level 
diagram by the multiplicative factor 
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V 



X 




Fig. 2 Tree level contribution to inclusive semi-leptonic B de- 
cay. 

where SZ;, = ^S|^^^, and — iE is the self-energy loop of a 6 quark with momentum p. 
This loop gives, for D = 4 — e dimensions, 



-iE 



3tt 



[(e - 2)|^Soo(p, 0, rrib) + (4 - e)mbBo{p, 0, m^) 



where 



and 



— 1 



J [(/ - A2 + ie][P - + ie] 

- + ln(47r) - 7 + 2 + in^ ^ 

e 



p°'Boo(p,0,mi,) = 



IJ,^ p^ IJ,^ 



—1 

7r2 J [{I - pY - A2 + ie] [/2 - m2 + ie 

4 



1 1 , , 1 ^ . rnr - p^ 

In 47r -7 - l-f m ^ / 

e 2 2 ' 2p4 



+ 



m4 m| 



2p4 



2p^ 



+ 



4 4 2 

p — m? , p 
In 



2p^ 



mr 



1^' 



(4.2) 



(4.3) 



We may set the gluon mass to zero for the self-energy, which is infrared finite, but 
for the derivative of the self-energy the gluon mass A is needed to regulate the infrared 
divergence. In this case, we get 



d_ 



a 
3^ 



In 47r + 7 - 4 4- In 



m 



\2 

o 21n^ 
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Fig. 3 Final quark self-energy diagram. 



The resulting form factors are 



J.0 



rpO 

-'-2 



2Ao 



{rrib - v-q)(l + 6Zb) 



^0 



rpO 

-'-5 



1 



2Ao 




1 



2A, 



(4.4) 



where Aq = {p — — m'j + ie. 



4-2. Final Quark Self Energy 

Since we are using the optical theorem, the light c or u quark self-energy is not 
incorporated into wave-function renormalization, but instead arises as a separate graph. 
These contributions to the form factors enter in the same proportions as the tree result, 
since self-energy corrections don't alter the Lorentz structure of the leading order operator. 
The self-energy graph of fig. 3, which also accounts for final state radiation, gives 



(Ao) 



[2m2(4-e)So(p-?,0,mj)-f 
((p _ qf + ^2)(e - 2)Soob - 9, 0, m,)] . 
11 



It gives a form factor contribution of 



-,self '^s 1 



where 

/ 2 \ 777,^ 

^[(p - m?] = - - 7 + ln(47r/i') (-(p - qf + 7m?) - (p - qf + 10m? ^ 



6 4 \ 

777/ ■ 777 ■ \ 

(^-3^ - - (P - + 7m? j ln((m, + A)^ - (p - - ie) 

(4.6) 

The optical theorem tells us that to obtain the B decay rate we need to compute the 
imaginary part of the T^'s. In the case of the heavy quark self-energy, as 5Zi, is real, there 
is only one contribution, coming from 

Im{^^} = -Tr5{x). (4.7) 

{x + le) 

On the other hand, Im{T^'^^^) has two contributions, a first one from 

M(^} = rfW, (4.8) 

and a second one from the imaginary part of the function A[{p — q)'^,m'j]. Those two 
contributions correspond to two different cuts in the Feynman diagram. In the first, we 
cut the final quark propagator near the weak vertex and put it on-shell with the Dirac 
delta. The virtual gluon renormalizes the light quark wave-function. It is easy to see, 
by integrating the S and S' functions over v-q, that the resulting initial and final quark 
wave- function renormalizations are identical, after the obvious interchange m^ <-> rrij. The 
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Fig. 4 Vertex correction diagrams. 

second contribution corresponds to cutting through the loop. In this case, the gluon is 
real, corresponding to final state bremsstrahlung. The fact that A only has a non-zero 
imaginary part when {p — q)"^ ^ ("^j + A)^ immediately implies the parton model branch 
cut ( |2.7| ). A has an additional dependence on the gluon mass which cancels the infrared 
divergences that arise from the vertex corrections. The analogous infrared divergence from 
bremsstrahlung off the initial quark is contained in the box diagram, fig. 5. 



4-3. Vertex Corrections 

The virtual gluon correction to the weak vertex is 



3 J (27r)^ [(/-p)2-A2 + ie][(/-Q)2-m2+ie][Z2-m2 + ie]' 
This enters into both vertex graphs in fig. 4, giving 



(4.9) 



where 



and 



SttAo 



^[.,a,af3] ^ 



7r2 J [/2 -ml+ ie] [(/ - - + ie] [(/ - - + ie] 



J 

2 

J 



'PlIx 



\'j()')YPL + '^^'^olIx'P'YPl 



N 



{V) 
af3 



= l^lcYllilxfYPL + l^i>'l^lcYPLll3lX, 
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(4.10) 



(4.11) 



(4.12) 



with p' = p — q. The Dirac algebra from eq. ( 4.12 ) is dramaticaUy simphfied by using 



heavy quark symmetry. Spinors Ub{p) correspond to operators b = h^^^ + C( — ,), so to 



leading order in As£iZ. have 

{B\u,ip)Tui,ip) \B) = ixrpi^]. (4.13) 

The same formula applies to B* mesons if we average over spin. Higher order corrections 
in can be computed by inserting heavy quark symmetry breaking terms from the 

heavy quark Lagrangian fl^ . 



The integrals C*! '"'"^! may be expressed, by means of the usual Feynman parametriza- 
tion, momentum shifting, and integration, in terms of the Feynman parameter integral 



5'(a, 6, c, (i, e, /) = dx dy In X{x,y; a), (4-14) 

Jo Jo 



X{xj y; a) = — {ax + by + cx^ + dy^ + exy + /), 



where 

1 

and 

a = -{p^ -ml + A^) 
b = — }? — {p — q)^ 
c = p^ 

(4.15) 

d= {p- qf 
e = 2p ■ {p — q) 

The real and imaginary parts of this function are evaluated analytically in Appendices 
A and B. The special case of a massless final state quark, which would give spurious 
divergences if the zero mass were substituted in the formulas in the previous appendices, 
is considered in Appendix C. These closed expressions are exact, and have been verified 
numerically. However, they are cumbersome to handle, and it may be simpler to evaluate 
these functions numerically than analytically. Expressions for the C integrals in terms of 
S and its derivatives are given in appendix D. 
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After evaluating the traces and integrals in terms of S and its derivatives, we arrive 
at the following results: 



rpvert 



Twert 
2 



rpvert 



^^[(-. + v.,)(i-| + iln4.-l-5) 

+ mbm''^Sj - ml i2S,a + 2S,b - S,^ - S,d - 2S,, - Sj) 
+ mlyqi3S,a + 55,6 - S,, - 3S,d - 4-5,e - 2Sj) 

+ nibq^ {S^a + S^d - Sj) + q^Y -q {S^b - S^d) 
- 2m6(vq)2 {S,a + 2S,b - S,d - S,e - Sj)] = 



a. 



SttAo 



vert 



rjivert 



«s 4 . ,1 7 1 , , , ^^ 

+ {ml + m^m] - m^f^) {S^a + S,b - Sj) + mbq^S^e)] 

cts 2 . , 1 7 1 , , , „v 
— — - - - + - ln47r - 1 - -S) 
37r Ao e 2 2 ' 

+ ml {-iS.a + 38 ^b - - S,d - 25,e - 2Sj) 

+ g2 (5 6 - S,d) - 2mbV-q(5a + 2S,b - S,d - S,e - Sj)] 



SttAo 



rpvert 
-^2 



a. 



SttAo 



rpvert 
^3 



-as 4 



rpvert 



-r— mb (S a — S e) 



-as 2 ,,1 7 1 



T^wert 
4 



SttAo 
1-5) 



— [(--- + -ln4TT 
Stt Ao ^^e 2 2 

+ (5 a + 25,6-25/) +m|5,6 
- 2m6V-q(5a + 5,6 - 5,e - 5,/)] = 



-n 



vert 



(4.16) 



SttAq 

The notation 5,a means the derivative of 5(a, 6, c, d, e, /) with respect to a. Again, 
we have two contributions to Im{T^^'^^): one from the imaginary part of the light quark 
propagator and one from the imaginary part of 5 and its derivatives. The former corre- 
sponds to the virtual gluon corrections to the vertex, and the latter to bremsstrahlung due 
to interference between initial and final state radiation. 



4.4- Box Diagram 

The box graph, depicted in fig. 5, gives 



d^; 



7-(/ + mbh^'PLi/ - + mbhc 



(2Tr)^ [(/ - p)2 - A2 + ie] [(/ - q^ - m] + ie] [/2 - + ie]2 
Stt 



«6(p) 



(4.17) 



Ub{p), 
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Fig. 5 Box diagram. 



where 



jj[.,a,al3,aj3^] _ 



ml + ie]2[(Z - g)2 - m2 + ie] [(/ - p)^ - A2 + ie] 



(4.18) 



and 



^if = rribl^-i^l^PLlpYlx 

+ rrib'-i^'-i^PLlcY'lp'lx - l^lcl^PhiYlplx 
= I^IcI^PlIpYi^Ix. 
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(4.19) 



Evaluating the traces and integrals as with the vertex diagrams yields 

OTT 

+ S^a,d + 2S^a,e + ^S^aJ + S^c,c + 3S',c,d + SS ^c,e + S ^d,e) 

+ mbq^{-2S^a,b + S^a,d - 2S^a,e + S^c,d + S^d,e) - V-q {S^a + 5',e) 

+ ml^-q {2S^a,a + QS^a,b + S^a,c - 5',a,d - 2S^a,f - 65',c,d - 35',c,e - ^S^d,e) 



+ q\-q {S,a,d - S,d,e) + mb{vqy{-2S,a,d + 2S,c,d + 2S,d,e)] = 

OTT 



+ mU4S^a,a + 4^,„,, - S,a,c ' ^,a,d - 2S^a,e ' 2S, 

— 3S',c,d — 3S',c,e — S^d,e) — 'mbq^{S^a,d — S^c,d — S^d,e)] = —^T^°^ 

OTT 



OTT 

Wj, ( 2S^a,a 2S^a,b S,a,c "I" S^a,d "I" 2S^a,f 2S^c,d ^,c,e S,d,e) 

- q\S,a,d - S^d.e) - m,vq {-2S,a,d + 25,,,d + 2^,^,,)] = ^f^^^ 

OTT 

^box ^ 4 [-mb{S,a,b + S,a,e ' ^,c,d - S ,d,e) + V-q {S,a,d - S,d,e)] ^ ^ 

OTT OTT 

OTT 

- ml{—2S^a,a — GS^a,b — 5',a,c + S^a,d + 2S^a,f + QS^c,d + 35',c,e + 35',d,e) 

- q\S,a,d - S,d.e) - mbyq{2S,a,d - 2S,c,d - 2S,d,e)] = ^f^"\ 

(4.20) 

The function S^a,b denotes the second derivative of S{a, 6, c, d, e, /) with respect to a 
and b. There is only one way to obtain an imaginary part from the T^°^ , via the imaginary 
part of 5" and its derivatives. This corresponds to cutting through the loop to obtain the 
bremsstrahlung due to initial state radiation. 



5. Analytic Structure of Form Factors 

The form factors VFi(v-q, g^) are given by the discontinuity of Ti(v-q, g^) across the 
cut in the complex v-q plane. The general location of the physical cuts may be deduced 
from the expression ( |1.6| ) by inserting a complete sets of states. There is a cut along the 
real axis for 

m R + q^ — m\ 
v-q<— ^ ^. 5.1 

^- 2mB 
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This cut corresponds to weak amplitudes for which the initial hadronic state is a S (con- 
taining one b quark and one light antiquark), and the final state contains one j type quark 
{u or c) and the same light antiquark. mx is the mass of the lightest such hadronic fi- 
nal state. The weak decay B XrVj- corresponds to the portion of this cut for which 
v-q > ^/q^ (since physical final state lepton pairs have > 0). By crossing symmetry 
one obtains cuts in other portions of the above region, corresponding to different physical 
processes, such as Bur Xt. 

There is one other class of cuts arising from the a;*^ < portion of the integral in 
equation (|1.6|), and corresponding to the production of two h quarks in the final state. 
This cut lies along the real axis for 

v-q> ^ \ (5.2) 

where m'x is the mass of the lightest hadronic state containing two b quarks which can be 
produced in a semileptonic weak process involving a i? in the initial state. This process 
occurs only at order a^. 

The location of the branch cuts for the form factors Ti in the complex v-q plane for 
fixed can be explicitely determined by looking at the discontinuities of the T's along the 
real v-q axis. For the vertex and box diagrams, this means looking for an imaginary part 
of 5'(a). Im [5'(a)] is simply the area within the triangular integration region over which 
the argument X(x,y;a) is negative, the branch cut is determined by studying the curve 
X(x,y;a) = 0. After expressing a in terms of physical variables, it is easy to show that 
Im [5'(a)] 7^ if and only if {p — q)'^ > m?-. This implies a branch cut for v-q < — — 
in agreement with ref . For the final quark self-energy, we have to look for an imaginary 
part in A \{p — qY, m'j~\ , which leads us to the same cut. 

An alternate way to see this is to consider two particle unitarity. Unitarity of the 
amplitude Aif between initial state \i) and final state (/| implies 

n 

For a two particle intermediate state (in our case a gluon and a quark), the right hand 
side becomes 

i k 

where k is the magnitude of the three-momentum of either intermediate particle in the 
center of mass frame, and Et is total center of mass energy. Since 



J dn (/I A \ki, k2) {ku k2\ A^ \i) e{El - mj). 




k I rn^j 



Et 2\j' {p-qY' ' 
two particle unitarity provides a positive check on our calculation. 
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6. Divergence Cancellation 

6.1. Ultraviolet Divergences 

We performed the loop momentum integration in 4 — e dimensions in order to regulate 
the ultraviolet divergences. The divergent part of each graph is always proportional to the 
factor 

2 

- -7 + ln(47r;u2). (6.1) 

e 

The final and initial quark self-energies contribute the same ultraviolet divergence, while 
the box is convergent. One can see from ( [4.4| ) and ( |4.16| ) that the divergent contribution 



from the vertex is -2 times that of the initial quark self-energy, so the dependence on the 
regulator e and the renormalization scale /i cancel exactly. This is nothing but an explicit 
application of the Ward identity to our case. 

6.2. Infrared Divergences 

The cancellation of infrared divergences (which we have regulated with a gluon mass) 
is more involved. The triple differential decay distribution is IR divergent; the dependence 
on the gluon mass cancels only after integration over v-q. The self-energies exhibit the IR 
divergence explicitely as In terms, but in the vertex and box diagrams the divergences are 
hidden in the S functions. It can be seen by numerical evaluation that of all the S functions, 
only Sj and S^aj blow up for zero gluon mass. These two functions can be written as 
a gluon independent term plus another term proportional to InA^. Upon integration 
over v-q, we have verified that the coefficients of these logarithms cancel exactly the IR 
logarithms from the self-energy. This non trivial cancellation of the gluon mass dependence 
give us confidence in our formulas. We have numerically verified this cancellation for up 
and charm final quark masses and for several values of . The cancellation is independent 
of the tau energy, since all the T form factors are independent of E^-. 



7. Ab baryon decay 

We now calculate the additional spin-dependent form factors Gi through Gg that 
describe the decay of a polarized A^. The Dirac algebra is a little more complicated in 
this case since we cannot use equation eq. ( |2.9p . First we reduce the expressions given in 
eqs. ( ^.121 ) and ( |4.19D by means of the Chisholm identity 

7 7 1=9 1+9 i^-g^ 1 +te^ ^Iplb- I'.lj 
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Note that in our convention £0123 = +1- After some algebra which we performed with 
the help of Feyncalc [|TB[ the Dirac structure of the N's reduces to terms with zero or one 



gamma matrix and 75 's. The matrix elements of these operators can be further simplified 
using the heavy quark symmetry. To first order in l/rrih we have 

{Kb{v, s)\hbThb \Ab{v, s)) = u{v, s)Tu{v, s), 
u{v, s)u{v, s) = 1, 
u{v,s)'j5u{v,s) =0, (7.2) 
u{v, s)'y'^u{v, s) = v'^, 
u{v, s)7'*'75w(f , s) = s'^. 

We now contract the resulting expressions with the C's given in Appendix D and collect 
terms with respect to their Lorentz structure. The baryon form factors Ti to T5, which 
enter into the first line of equation ( |3.9| ), are identical to the meson form factors Ti to T5 
in equations (|4.4| ), ( ^75|) , (|4.16| ), and (|4.20[) . The spin-dependent form factors in the rest 



of (|3.9| ) are found by comparing the above result to the defining equation for the G's, eq. 
We divide the spin dependent form factors into three pieces 



Q ^0 _|_ Qvert _|_ Qhox 

The results from wavefunction renormalization and the self-energy graph are proportional 
to the tree level result, so they may be combined into 



-1 

-mi 

2Ao 
1 



2Ar 



(1 + 5Zb + 






m] 


'-{l + 5Zb + 






m 


(1 + 5Zb + 








(1 + 5Zb + 






m] 


(1 + 5Zb + 









2Ao' SttAo ^' (7.3) 

mb 

2Ao 
-1 



2 — (jrg — (jr4 — (jrg — U, 



where the functions and A are defined in eq.s ( [4.4| ) and ( |4.6| ). 
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The contributions to the G"s from the vertex diagrams are 



" 2 + ^ ^''^^''^^ + 1 + ^ - ml{3S,a + 35,6 - - S,a 
- 2S,, - 2Sj) - q^{S^b - S,a) + 2m,vq(5,„ + 25,6 - 5,<, - 5, ^ - 5, j)) 

revert _ '-^ Q 
— 1 — "J e 

^ SttAo ' 



Invert 
<J"4 



=0 



revert ^ 2Tn{, , , 



3n Ao 

=^^^-^1 ^^^^^^^ + 1 + S- m'^{S,a + - Sj) 

- ml{S^a + S^b - Sj) + q^{S^a + S^b - - Sj) 

- mby-q{S^a - S^c - ^e)) 
^vert ^«2^1_2^1 ^^^^^^ _i_s + m]S,b- m,vq(5,„ + 25,6 - 5, e - 25, 

+ m^(25,„ + 25,6 - 5,c - 5,e - 25,/)) 

vert (X rtlb ,^,\ 7,1 



GT' =T^^(2(- - + § ln(47r)) - 2 - 25 - q\2S^a + 25 ^ + 5 e - 25,/) 
OTT Aq e z z 

- 2m25,/ + 2m^(25,„ + 25,6 - 5,^ - 5,^ - 25,e - 5,/) 

- m6V-q(35,a + 65,6 - 5,c - 45,^ - 55,e - 25,/)) 
Qvert 1 /_2(1 -2 + 1 ln(47r)) + 2 + 25 - 2q^(S 6 - 5d) 



- m2(35,„ + 45,6 - 5, - 25,^ - 35,^ - 25 /) 
+ m6v-q(45 „ + 85,6 - 45,^ - 35 ^ - 45,/)). 



(7.4) 



The contribution to the G's from the box are calculated in the same way as the vertex 



contributions. Again, the algebra is a little tedious but straightforward. Making repeated 



use of the Chisholm identity eq. (|7.1| ) , the heavy quark matrix elements eq. (|7.2| ), and 
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contracting with the -D's eq. ( 4.18 ) we arrive at 




'box 



a 



a 



box 



a 



^5°^ =7^2m6( — 5" a,& + 3S a,d + S a,e — '2S c,d — '^S d,e) 

on 



OTT 

+ 3S'^c,e + 'S'^d.e) + '^b(^{^S ^a,b — 4:3 ^a,d — S^a,e + S ^c,d + S ^d,e) 

- SrUbiS^a - S^c - S^e) - 'mlv-q{3S^a,a + 7S^a,b - S^a,c " 7S^a,d " 8S^a,e 

— '^S^aJ + 4S^c,d + 25',c,e + 25',d,e)) 

^9"^ —-^iS,a — 3S',e + q^{S^a,d — S^d,e) — ■q{S^a,b + S^a,d + S^a,e — '^S^c,d 



The self-energy diagram gives no spin dependent contribution. 
8. Conclusions 

We have presented as corrections to triple differential inclusive distributions for B 
XqTV, in closed analytic form. We have also calculated the dependence of the distribution 
on the tau polarization, for both up and charm quarks. We have not performed the 
integration over neutrino energy which is necessary to give the physically meaningful double 
differential distribution, but have instead left the answer in terms of a one dimensional 




2S^d,e) + ml{S^a,a + S^a,b 




(7.5) 
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integral which can be evaluated numerically. An additional integration over tau energy or 
momentum transfer will yield ^ or -^-^ respectively. Formulas appropriate for insertion 
in a Monte Carlo program are presented in Appendix E. 

We have verified numerically the cancellation of infrared divergences due to soft gluons, 
and the stability of the double differential distribution against variations in the infared 
regulating gluon mass. In addition, the formulas for the graphs and integrals entering the 
triple differential distribution have been checked extensively numerically. 

Our computation has been in the context of a heavy quark operator product expan- 
sion. An advantage of this is that extending the calculation to polarized A's is simple. 
The quartic distribution for polarized Lambda's to polarized tau's +z7 + Xq involves nine 
additional form factors, which we have computed. Again, integration over the neutrino 
energy is necessary to get an experimentally observable distribution. 

The formulas in this paper may be useful for purely theoretical reasons, in addition to 
experimental ones. A computation of the corrections to the Wilson coefficeints of the 
operator product expansion may give insight on how perturbation theory breaks down as 
one nears the endpoint of the lepton spectrum. This requires calculating one-loop graphs 
in QCD (which we have done) and comparing them to one-loop graphs in the effective 
theory. 

Strong corrections to form factors at 0{-^) may also be of interest [0. The structure 



of our results is amenable to such a calculation, by expressing the h momentum as p = 
-|- k and expanding in residual momentum k. The only complication here is that care 
must be taken not to cut the b quark lines when computing Bremstrahlung; this requires 
a modification of our formula for the imaginary part of S. 
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Appendix A: Evaluation of Re(S) 

We want to evaluate 

S{a,h,c,d,e,f)= dx dy In X (x , y; a) 

Jo Jo 
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where 

X{x,y; a) = —{ax + by + cx"^ + dy^ + exy + /). 

We evaluate the real and imaginary parts of S separately because we will use a change of 
variables which causes our integration contour to cross branch cuts for some values of our 
parameters, producing spurious contributions of m. In addition, this change of variables 
is valid only ii d ^ For an on-shell up quark, we may either take the limit as (i ^ 0, or 
use the formula in appendix C. 

First, we perform a change of variables, x ^ 1 — x and y ^ y — Oix, where 

a= -^{e+ \J — 4cd), 

and use the identity 

"(1— a)a; i-l—a nl 



pi p{l—a)x pi— a pi p—a. pi 

dx dy — dy dx — dy 

Jo J-ax Jo J-rr^ 

(1 — a) ot 



dx. 



In the first integral on the right hand side, make the replacement y — > (1 — a)y, and in the 
second, replace y — — ^- The resulting equality is valid only for the real parts: 

Re[S{a.)] = Re{ [ dy f da; [(1 - a) InX(ai) + alnX(a2)]}, 

Jo Jy 

where 

— {(^ij ^i-i (^ij Cj, /j}, 

ai = a{b + e) — a — 2c = 02, 

61 = (1 - q;)(6 + e), b2 = —a{b + e), 

Cl = = C2, (^.1) 

di = {1 — a)^d, d2 = o?d, 

Cl = (1 — a){2da — e), 62 = a.{e — 2da), 

fi = f + a + c = f2. 

This expression can be evaluated straightforwardly by completing the squares inside each 
integral. The final result is 

Re[S(aL)] = Re[(l - a)F{aLi) + aF(aL2)] iA.2) 

where 

F(a) = - [-/(as) + /(a4)] , 
e 
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/(a) = a{g{a, b, c, d) + g{a, —b, c, d) — 1 + 3c — 2d 

+ - c + In a + (c^ - - In a) [ln(c + d) - ln(l + c + cZ)]}, 



(A.3) 



5f(a, 6, c, d) = ln(6 + d) 
1 1 



I52 + 5c + cd - ^ ^^2 _ in(^) 

^ ^ CO 



+ ln(l + 6 + d) 



^ - - c - 6c + d - cd + + (c^ - 6^) in(l±^) 
2 2 2 c — 



(^.4) 



03 = d + e, a4 — d, 



f 



2d + 2e^ d + e' 



b+e a+f 



2d 



d 



a a + b 
C3 = - 



a 6 + e 

e 2d + 2e' " e ~ 2d ' 



d; 



a + 6 



d4 = 



6 + e 



^ 2d+2e' 2d 
63 = 64 = /s = /4 = 0, 



and the Spence function is defined so that 



Sp{x) 



° ln(l - t) 



dt. 



(A.5) 



Note that in the evaluation of, for example, F(ai), the parameters which enter into the 
definitions of a3 and a4 in equation (A. 5) are ai. 

Appendix B: Evaluation of Im(S) 

It is much easier to evaluate the imaginary part of S than S itself. 

ImS = dx dy ImlnX{x,y;a) 

Jo Jo 

is just TT times the area in which X is negative, inside the triangle in the {x, y) plane defined 
by the integration limits (recall that / = — ie is the only complex input parameter). 
To obtain this area, we solve the equation X = 0, which gives us a hyperbola in the (x, y) 
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plane, and integrate x as a function of y (or viceversa) in the triangle. The evaluation of 
this area gives: 

ImS = 7r(^(Vfe2 _4^y) _ A{-y^^~Idf)), 

where the function A is 

a{h~K) e{h-Kf f 2 6c- ae -h + K 



^^^^ Acd 16cd2 V-4c(e2-4cd) ^ 8cd )^ 

2a'^d'^ -2abde + b'^e^ -2de^ f + 2adeK -be^K 

{abe-b^c-a^d- f{e^ -Acd)) ^ _(^e^ _ 4cd) {b - K) 



(e2-4cd)^ d 

r-^ — - Ua'^d'^ + 2(2ade-be'^)(K -b) -Ade"^ f, 

Ve^-Acd^ ^ ^ 

(5.1) 

All the derivatives of S needed for the decay rate can be obtained from this formula. 

Appendix C: Evaluation of Re(S) at threshold, for massless final quarks 

The formula for the real part of S breaks down for d = 0, and it is cumbersome to 
take the limit numerically. It is convenient to have available a simpler representation of 
S which is valid when d = 0. The evaluation is straightforward, so we omit details and 
present only the final result: 

Re[S{a, 6, c, 0, e, /)] = Re[~ + ^ ln(c - e) 

- I ln(c)[i + 2ai - ^ + (c? + {-a, + V) Hj;^)] 
el e e b + e 

+ - ln(c - e)[^ + a2C + ts + U ln(i±^^)] (C.l) 

+ -/o(a2, — , C2) - -/o(ai, oi - -, ci) - /o(a2, -1 - 02, C2) 
e c e e 

+ ^i5'o(ai, -ax + -, ci) + —go{a2, ^2, C2)] 
e e 



where 



/o(a, 6, c) = /oo(a, 6, c) + /oo(a, b, -c) - ^ - a - 2b, 
go{a, b, c) = goo{a, b, c) + goo{a, b, -c), 
/oo(a, b, c) = (-^a^ -ab-bc+ ^c^) ln(a + c) 

+ {^ + a+^a'^ + b + ab + bc- ^c^) ln(l + a + c), 
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(C.2) 



and 



5ioo(a, 0, cj = in( — ^ — - — J in(l + a + c) - in(^ — -j in(a + cj 
a 

a — 6 + e 

ci 



(C.3) 



2(c-e) 



4c2 c 



(a-6 + e)2 6 + / 
"^ = ^/ 4(c-e)2 (^-4) 

^i = -[c?-(--ain 
e e 

^2 = 02 

e 

^3 = a — 2ca2 — ct2 

^4 = + (2ca2 — a)t2 + f — + 002- 

The real part is taken because, as in Appendix A, the changes of variable contribute 
extraneous factors of Itt to the imaginary part of S. While both equations (A. 2) and (C.l) 
are exact, single precision numeric evaluations of them may work poorly when the inter- 
mediate parameters are anomalously large or small (for example, when c ~ e above). 

Appendix D: The loop integrals in terms of S's 

C =-Sj 



+ (2 S,a + 2 5 6 - S,c -S,d-2 S,e - Sj) / 
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D =-S,a,f 



—S,a,b + {S,a,a + S,a,b ' S,aj) 

-{S,a,b-S,a,d-S,a,e) {q^ + q'' ) 

+ (2 S,a,a + 2 S,a,b - S,a,c - S,a,d - 2 S,a,e - S,a,f) / 

jja(3y = - 1 5 e (^^^ + ^"^^ + g""^ q^ ) - S,d,e q" q^ q^ 

- ^ iS,a - S,, - 5 e) {g^^p^ + g^^p'^ + ) 

- {S,a,d - S,c,d - S,d,e) ( g" g^ / + g^ g"^ + g" g^ ) 

~ {S,a,b ~ 2 S'^a^d ~ 2 ^a^e + 2 S^c,d + 'S',c,e + ^,d,e) 
( p'^ + qf^ p"" p^ + 9^ p'' p^ ) 

+ (3 S^a,a "l~ 3 5*^0,6 ^ 3 S^a,c 3 S^a,d 6 S^a,e 

- S,aJ + S,c,c + 3 5e,d + 35e,e + ^,d,e)p"/p' 



Appendix E: Double differential B meson Decay rate 

We collect here the complete formula for the ag correction to the double differential B 
meson decay rate. It is expressed as a one dimensional integral amenable to insertion into 
a numeric integration program incorporating any relevant cuts and detector efficiencies. 
The masses mr,mb,mj correspond to the tau lepton, the b quark, and the w or c quark 
masses, respectively. Ej. and Ejy are the tau and neutrino energies, v is the velocity of the 
B meson, and v-q = Er + E,^ in the center of mass frame. 
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The double differential rate is 



|2 ^2 



dq^dEr 

(i v-q( 



{nil, — v-q)-B[v-q, m 



(q2 _ ml) ((m, - vq,J(-4 + 3 Inm^ - 2 In A^) + fr%h 



+ 



2{ml + q"^ — 2mbV-q — m^)^ (m^ + q'^ — 2mbV-q 



+ 



n 



+ 



rfv-q (^2£;^(vq- E^) - 



2mb(-4 + 31nm^ - 2 In A^) + Ta"^''* 
2^ ' mf, i3[v-q, m^] 



{ml + — 2m6V-q — m^)^ 



Im{T^°^) 



+ 



{ml + q'^ — 2mbV -q — m 
+ ((2E, - v-q|t;,)(Q2 - ml) - 2ml{vq\th - E^)) ((-4 + 31nm^ - 21nA2) + 



vert I 



\th 



dvq{{2E^-vq){q^ - 

/m(t3^^'^*)/7r 
{ml + q'^ — 2m6V-q — m^, 

( q^ -ml \ 2 



m^ 



2ml{Y-q-E,)) ( 



i3[v-q, 



2(m^ + — 2mbV-q — m 



2^2 



+ 



+ 



Im{T{ 



TT 



-!vert 
-4 



(iv-q(- 



Im{fl^''^)/Ti 



ml + q^ — 2mbV-q — m^^ 



+ 



Im{fl'' 



TV 



+ 2ml ( (v-q|t/. - ^.)(-(-4 + 3 Inm^ - 2 In A^) + Tr%h) 

— i3[v-q, m^] 



- / d\-q{Y-q- E^) ( 

Im{f^^'^^) /tt 



+ 



(m^ + — 2m5V-q — m 



2{ml + q'^ — 2mi,v-q — m^)^ 

. ')) 



+ 



Im{fi 



'box\ 



TT 



where 
;B[v-q, 



q +2mbV-q + 7m- 



(m^ + — 2m;,v-q)2 {ml + — 2mbV-q) 
The T|t/i functions are the real parts of the T functions defined in ( [4.16| ) and ([4.20|) , 



evaluated at the parton model theshold v-q| 



th 



{ml + q'^ — m^)/2mf,. The real and 



imaginary parts of T may be expressed in terms of the real part of the function S evaluated 
in Appendices A and C, and the imaginary part of 5" evaluated in Appendix B. Integration 
limits are given in (3.5)-(3.7). Note that the T functions have an imaginary part only for 
v-q < {ml + q'^ — {mj + A)^)/(2mb); this A dependence of the integration limits cancels all 
gluon mass dependence upon integration over v-q, in the limit A ^ 0. 
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